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Abstract 

We compute correlation functions of inverse powers and ratios of characteristic polynomials 
for random matrix models with complex eigenvalues. Compact expressions are given in terms of 
orthogonal polynomials in the complex plane as well as their Cauchy transforms, generalizing pre- 
vious expressions for real eigenvalues. We restrict ourselves to ratios of characteristic polynomials 
over their complex conjugate. 

PACS: 02.10.Yn = Matrix theory 
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1 Main Results 



The theory of random matrices has found many applications in different branches of physics as well 
as in mathematics. One possibility to study the correlation functions of matrix eigenvalues, which can 
then be mapped to various physical or mathematical quantities, is to compute ratios of characteristic 
polynomials as their generating functional. Such correlation functions can also be studied in their 
own right, as they enjoy a direct physical interpretation as well. While the most general generating 
functional is known for real eigenvalues much less was known until recently for complex eigenvalues. 
Characteristic polynomials of the corresponding non-hermitian operators play an important role for 
example in scattering in Quantum Chaos, as reviewed in 2 , or in Quantum Chromodynamics 013]. 

Our purpose is to generalize the result for arbitrary products of characteristic polynomials 
and their complex conjugates of not necessarily the same order, to ratios of such objects. We restrict 
ourselves to the case of matrix models with a complex eigenvalue representation, to make the technique 
of orthogonal polynomials available. Its difficulty in the complex plane is that in general neither the 
three-step recursion relation nor the Christoffel-Darboux formula hold in general. We will still be able 
to show that part of the results of [H] generalize to complex eigenvalues, where streamline proofs of 
the previous achievements [7J |H1 for real eigenvalues are given. 

The partition function of a complex (matrix) eigenvalue model is defined as 



where we have introduced the Vandermonde determinant A^({z}) stemming from the Jacobian of 
the diagonalization. We suppose that the probability measure dw(z,z) = dzdz w(z, z) can be written 
in terms of eigenvalues and factorizes, and thus that eq. is invariant under permutations of the 

eigenvalues Zi=\ jsj. The weight function w(z,z) depending both on z and z shall be strictly positive 
on the domain of integration D. D can be either a compact set or the full complex plane. Furthermore 
we require that all moments exist, f D dw(z, z)z k < oo for all k = 0, 1, . . . . Examples for such general 
weight functions and domains D are given in jSJ. Under such conditions a unique set of orthogonal 
polynomials in the complex plane can be introduced using the Gram-Schmidt procedure. The monic 
polynomials iTk(z) = z k + . . . of degree k follow from 



with their squared norms > being strictly positive. We note that in contrast to the weight the 
polynomials tt^(z) only depend on z and not its complex conjugate z. The corresponding Cauchy 
transform in the complex plane is given by 



Note that the pole e is an integrable singularity in the complex plane. If we want to allow for the 
hermitian limit to be taken we have to require e ^ D. In that case D has to become compact at 
least in the large- N limit. The Cauchy transforms h n (e) only depend on e and not on e. Another 
possible Cauchy transform, dividing n n (z) by z — e, is not needed as it can be expressed through 
the 7Tfc< n (z) and ho(e). Expectation values of observables of symmetric functions of the eigenvalues 
In = /n(zi, • • • , zn) can be defined as 
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Our objects of interest are ratios of characteristic polynomials Dn[/j] and their conjugates, 
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We can now state our main result, generalizing [S] (see also theorem 2.13 of [Hj). 



Theorem: Let {/tj}j=i,...,£ an d {^k}k=i,...,M be pairwise non-degenerate complex variables. For 
< M < N it holds: 
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(1.6) 



The following two special cases are worth to be mentioned. For L = and M/fl, that is for inverse 
powers only, we obtain a determinant composed purely of Cauchy transforms eq. (|1.3|) . generalizing 
the results of |Hj (see also theorem 2.10 in jSJ). In the opposite case, for L / and M = with only 
products, we partially recover the result of [3] in terms of polynomials only. The theorem as well as 
the special cases trivially carry over to the complex conjugate expressions. The limit of coinciding 
variables, e.g. ^ = fij, can be easily taken, leading to derivatives of the polynomials and Cauchy 
transforms. 



2 Proofs 

The proof will very closely follow the steps taken in [Hj. Due to the Heine- formula it is well known 
that orthogonal polynomials with respect to a given weight can be expressed through characteristic 
polynomials, 

(D N [fj]) w = 7nv(/i) • (2-1) 
In the following it will be useful to consider a generalized measure, 

dw^ (z, z) = ^=^_ 3 ~ Z } dw(z, z), £,m>0, (2.2) 

as well as the corresponding quantities eqs. (|l.ljl - (jl.41) 1 . The expectation value in the theorem 
is then proportional to "orthogonal polynomials" with respect to eq. (|2.2|l . 7r|^ 1 (/xl)- We will 
explicitly construct such polynomials by requiring 

J dw le ' m] (z,z) 7rf m] {z) z k = 0, j > A: > . (2.3) 

They can be interpreted as bi-orthogonal polynomials [S]. Let us stress however, that our result eq. 
(II. 6j) can be entirely formulated in terms of truly orthogonal polynomials and their Cauchy transforms, 
which form a bona fide scalar product. 

1 The superscript [0,0] is sometimes dropped, corresponding to the previous definitions. 
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Our proof goes in four steps. We consecutively construct the polynomials TTn'°\z), 7rn' m \z), the 
Cauchy transform of the latter hn {z), and finally Tin' m \z). In this way we first show the special 
cases M = and L = respectively, before arriving at eq. (|1.6j) . 

Step 1. Let us define for £ > 1 



n n (pe) ... n n+ e(m) 

7T n (z) ... 7T n+ i(z) 



for which it holds 



Because of <Zn'°'(/^j) = 0, j = 1, - • • ,£, the ratio r^-j|p4^^ is a polynomial of degree 2 n. It can thus 
be written as a linear combination of the polynomials 7Tq ...,n( z ) forming a complete set. Consequently 



/ dw(z, z) q [ ^ 0] (z) z j = , < j < n - 1 . 



(2.4) 
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(m-z)... {fit - z) 



z J = 0, for < j < n 
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In order to achieve a monic normalization we can expand eq. Q2.4JI with respect to the last row, and 
take z — > co to read off the generalized Christoffel formula in the complex plane 
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. (2.7) 



These polynomials were previously computed in [S]. The denominator is non- vanishing due to the non- 
degeneracy of the fj,j. Due to the relation eq. (|2.1j) for the general weight ir^ \[ij + i) = (-DjvlMj+lD^y.o] , 
and the observation that the expectation value can be written as a telescope product, 
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we can deduce eq. (|1.6|) for M = upon using eq. I|2.7|) . 
Step 2. Next we define 

^n— ro(^l) 



,[0,m] 
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which automatically implies 



dw(z, z) 
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(2.9) 
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The fact that the degree is n and not less can be shown by induction. 
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For < j < n we can decompose 



where p{z) is a polynomial of degree < n — m. Consequently 

[0,m] 



(2.11) 
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dw [0,m] (Z; g [0,m]( 2 ) ^ = y ak dw(z,z)^ ^ + / dw(z,z) q^ m] {z) p(z) = (2.12) 

Jd e k - Z Jd 



vanishes due to eq. (|2.1U|) in the first sum and orthogonality in the second term. In monic normalization 
we thus have the generalized Uvarov formula 
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(2.13) 



Step 3. Let < m < n. The Cauchy transform of eq. (|2.13|) for the measure dw^ ' m \z) can be 
expressed in terms of the h n (z) by writing 
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Only the term in e = e m +i is non-vanishing, and thus we obtain from eq. (|2.13|) 
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. (2.15) 



These expressions can be used in the following identity for an inverse characteristic polynomial 
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After decomposing the inverse product we have used the identity = |Aat_i({z})| 2 \\ (zn—z^) 

11 " A Zk k<N 



k<N 



as well as the permutation symmetry of the integrand to deduce this generalized Heine formula for 
the Cauchy transform. In order to apply eq. (|2,15|) we rewrite identically 
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From eq. (|2.16|) valid for the general weight eq. (|2.2|) we can conclude 
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writing the expectation value as a ratio of partition functions. It is a well known fact that the partition 



function eq. <fl.l|> can be expressed in terms of the norms eq. (|1.2jl . = iV! Y\j=o r v or equivalently 
[N — k) = Zjsr-k/{, r N-k^N-k-i)- Replacing this factor in eq. 1)2.18)1 and inserting it into eq. (|2.17j) 
all partition functions cancel and we obtain 
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Together with eq. (|2.15|) this leads to the theorem eq. 1)1.6)) in the special case of L = 0. 

Step 4. We can now give the polynomials with respect to the most general weight eq. (|2.2j) 
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holding for < m < n. If we define by q n ' (z) the determinant in the numerator it holds 
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This can be seen following the same lines as in the previous steps, and thus that eq. Q2.2U|) is correct 
including its normalization. In order to prove eq. 1)1.6)1 we decompose 
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Inserting eq. 1)2.20)1 and the previous result eq. (|2.19)1 from step 3 the theorem eq. 1)1.6)1 follows. 



3 Conclusions 

We have computed the correlation functions of arbitrary products of characteristic polynomials over 
arbitrary products of complex conjugate characteristic polynomials for random matrix models with 
complex eigenvalues. This extends previous results for only products of mixed characteristic polyno- 
mials and their conjugates jS] (see also |H|). From the result jS] we expect that more general correlation 
functions of mixed ratios will contain both polynomials and Cauchy transforms as well as the various 
kernels constructed out of them, as introduced in 8 in the real case. This would be needed to compare 
for example to the matrix model result (without eigenvalue representation) ^U], where the average of 
the inverse of a single characteristic polynomial and its complex conjugate was computed. We hope 
that our results will help to further clarify the issue of complex matrix model universality jllj . with 
an extension of the results [HIE] being desirable. 
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When writing up our results the preprint JB] appeared, which partly overlaps. There, ratios of 
mixed characteristic polynomials and their conjugates are expressed in terms of three different kernels, 
providing formulas of the two-point functions type. 
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